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1 . Introduction 


Due to screening by dust particles, the UV radiation field of 
inters teellar origin is practically inexistent within very dense 
interstellar clouds. However, as argued first by Prasad and Tar- 
afdar (1983), it appears possible that the cosmic-ray excitation 
of^ the Lyman and Werner systems of the hydrogen molecule could 
°^^-Si na Le a chemically- s ignificant flux of UV photons even within 
such dense clouds. 

This suggestion by Prasad and Tarafdar was investigated later 
by Lepp et al (1986), Sternberg et al (1987), Gredel et al (1987), 

found that cosmic -ray induced photon production may signifi- 
cantly affect the chemistry of dense clouds. 

In view of the important implications that the existence of 
such photons may have for the gas and dust evolution, and taking 
into account that the transfer of radiation inside a cloud is 
strongly affected by extinction and scattering properties of dust 
particles, as well as by their distribution, it seemed worthwhile 
to re-analyze the c.r.- induced production of UV photons, giving 
particular emphasis to the treatment of the transfer of the 
emitted photons. 

Elsewhere (Aiello et al, 1987) we computed the c.r. 
induced UV flux in a homogeneous spherical cloud as a function of 
the extinction and scattering properties of dust. In the present 
work we extended the analysis to the case of clouds with radial 
gas and dust density gradients. 

2. UV Photons production and transfer 

The equation of statistical equilibrium for the vibrational 
levels of the upper electronic states (B^u*, C^Ilu) is- 

14 14 

Huff j (A VXj + B uj <I uj >) = n o( C ou + B ou <I ou» + ^j n j B ju <I uj > [1] 

where n u , n Q , nj are, respectively, the populations of the 
upper electronic vibrational level, of the ground electronic level 
with v''= 0 and of the j th vibrational level of the electronic 
ground. 

Because of the low level of expected UV fluxes we can disre- 
gard the stimulated emission with respect to the spontaneous emis- 
sion. We can also overlook the absorption from electronic ground 
states with v''>0. 

Finally, since the photons emitted in the (u v' '=0) transi- 
tions are absorbed very close to the point at which they are gen- 
erated ("on the spot"), we can put the corresponding fields as 
equal to the source functions. 

So, the equation [1] becomes: 

1 4 

n u — n o^ou j J \ij ) 1 [2] 
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The transfer of emitted photons has been computed by using a 
method the outlines of which have been given by Van de Hulst and 
Davis (1961), adapted to a spherical symmetry with a radial gra- 
dient of dust density. 

The approach is the following: let I Q be the intensity of the 
zero order and let I n , the intensity of the light which has 
been scattered n times in succession. We will proceed step by step 
from the I n to the I n+1 , and will then compute the total scat- 
tered light: <I> = <I Q > + E n <I n >. The intensities are averaged 
over the line profiles as well as over the solid angle. 

The source function for the scattering of order n in each 
point of the cloud is related to the intensity of the order n-1. 
In other words the calculation scheme is: 

S Q (molecular emission) => I Q => => I L =* S 2 => I 2 ^ 

Mathematical and computational details can be found in Aiello 
et al (1984). 

For the photon fluxes we obtain the following expression: 

1 4 

N uj - c ou R F(r/R,i/ Uj ; r >7 ,g) Vj n D (l-tf) (ZA^ ) ' 1 [3] 

Let us to discuss the quantities involved in [3]: 

is a correction factor that allows for the fraction of 
absorption events that lead to the dissociation of the molecules. 
Following Jura (1974) we put \}> = 0.1. 

The excitation coefficients, C QU were computed by making use 
of the results from analysis of the effects of the energy deposi- 
tion by c.r. within a molecular hydrogen cloud performed by Aiello 
& Cecchi-Pestellini (1987). In this computation the interstellar 
proton spectrum by Webber & Yushak (1983), extrapolated down to 1 
MeV, was adopted. A check on the validity of our choice is repre- 
sented by the value of the resulting ionization rate (*=3x10" 17 
s“* per H 2 molecule, wich is in good agreement with the values 
derived from the analysis of thermal balance within dark clouds 
(Duley & Williams, 1984). 

The efficiency factor F is a function of the cloud structure 
as well as of the extinction and scattering properties of the 
dust . 

a) Extinction properties 

The wavelength dependence of the optical depth depends on the 
extinction law adopted. The latter has been studied in the direc- 
tion of many thousands of stars located in different galactic 
regions and as trophysical environments. The mean extinction curve 
(MIEC) so derived (Savage 6c Mathis, 1979) can be considered as 
representative of the general properties of dust in the diffuse 
medium and in old associations (Aiello et al . , 1988). However, 
extinction and polarization data show that in dense clouds and in 
regions of recent star formation, dust grains are likely to have 
different properties: in particular, their size distribution 
appears to be biased towards larger radii that results in a more 
or less marked flattening of the extinction curve in the far UV 
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spectral region. Therefore, the use of MIEC in radiation transfer 
computations could be completely inappropriate in the case of 
dense clouds. Because of this, we shall adopt here the extinction 
law derived for a-Sco (Snow & York, 1975), assumed as representa- 
tive of the extinction curve in a class of dense clouds. It is 
worthwhile noting that the a-Sco extinction curve is a conserva- 
tive choice. Indeed, the mean extinction curve towards inner 
region of p-Oph is even lower. 

b) Scattering properties 

7 and g can be derived, in principle, from observations of 
reflection nebulae and of the diffuse galactic light. However this 
task is more difficult to accomplish than the observation and 
interpretation of the interstellar extinction law. in fact only 
under exceptional conditions is it possible to derive the albedo 
and the asymmetry factor independently of one another. Further- 
more, the determination of 7 and g can be affected by the uncer- 
tainty of the geometrical relationship of the sources of illumina- 
tion and the scatters, as well as and by the uncertainty of the 
illuminating spectrum. The observational results are particularly 
ambiguous in the far UV region. 

Nor does the situation look much better when we turn to the 
existing models of the interstellar dust. Indeed, since the inter- 
stellar extinction curve is quite insensitive to the exact compo- 
sition of the dust, it can be matched by models which exclude each 
other. So existing models of dust based on the fitting of the 
interstellar extinction curve give different values of the scat- 
tering parameters for the same spectral interval. Moreover, as 
stated above, the properties of dust inside the dense clouds may 
differ from the properties of the dust in the general interstellar 
medium, so the uncertainty is even greater. 

It is possible, however, to derive from the observational 
results some indications about the scattering properties of inter- 
stellar dust. Indeed it appears that the dust is basically forward 
scattering and that the its albedo is smaller towards shorter 
wavelengths than in the visual (De Boer, 1986). Moreover, it 
appear possible, independent of particular grain models, to put 
constraints on the values allowed for the mean albedo and asym- 
metry factor of the dust responsible for the MIEC (Chlewicki and 
Greenberg, 1984). So it is possible to exclude from consideration 
the cases of very low values of the asymmetry factors coupled with 
large values of the albedo or of very large albedo coupled with 
completely forward scattering. Both the Mathis, Rumpl and Nordsieck 
(1977) and Greenberg and Chlewicki (1983) models give allowed val- 
ues . 

The validity of this analysis is limited to the grains as 
defined by the mean extinction curve. No scattering models, in 
our knowledge, have been worked out for the dust within the dense 
clouds. Some indications, however, can be obtained if the defi- 
ciency of small grains in size distribution is taken into account. 

In the present computations, we used the scattering parame- 
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ters measured in Orion by Mathis et al (1981) . The adopted data 
are those of their model 5. The reason of this choice is that 
Orion, being a region particularly suitable for investigating the 
scattering properties of the dust, exhibits in the same time a 
flat extinction. 

Finally, j , the probability for the 629 transitions consid- 
ered in the present work are from Allison and Dalgarno (1970) . 

3. Results 

We carried out computations of photon fluxes for two differ- 
ent model of radial density distribution inside the cloud (gas and 
dust densities are assumed to have the same density distribution) : 

I. The cloud has a small uniform density core (r c - 0.1 pc, 
n(H 2 ) = 10 4 cm“^) plus a r“l density gradient. The resulting total 
visual optical depth (edge to center) is 10. 

II. The cloud has a small uniform density core (same parame- 
ters as in the case I) plus a r“^ density gradient. The total 
visual optical depth is 6. 

In both cases the cloud radius is assumed to be equal to 1 

pc. 

The assumption that the interstellar UV radiation cannot 
penetrate inside the cloud obviously implies the existence, around 
our object, of an extended envelope, or that the object is located 
in a less dense but larger region (an example of such objects is 
the dark cloud B3 (Stenholm, 1985)). 

The photon fluxes obtained at the center (C) , in the middle 
(M) and at the edge (E) of the cloud are given in table 1. The 
resulting line spectrum is shown in fig. 1. 


Synthetic Spectrum 


center, case i 



Fig.l Synthetic emission spectrum in the Lyman and Werner bands 
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Our results strongly support the suggestion that the low 
energy cosmic rays may generate, inside dense interstellar clouds, 
UV radiation fields that may significantly contribute to the gas 
and dust evolution. 


Table 

i. 

PHOTON 

FLUXES (cm' 2 s' 1 ) 

x 10 

- 1 




CASE 

I 


CASE II 



c 

M 

E 

C 

M 

E 

Lyman 

275 

220 

115 

470 

225 

80 

Werner 

300 

240 

125 

515 

260 

90 
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